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In the minimal SU(3)l x U(1)n gauge model with a global L e — L M — L T (= L') symmetry and a 
discrete Z4 symmetry, it is found that the interplay between neutrinos and charged leptons contained 
in triplets of ^=(14,, t L , if) (i=l,2,3) naturally leads to the large mixing angle (LMA) MSW 
solution. The model includes two (anti)sextet Higgs scalars, with L'=0 and <S , ' + - ) with L'=2, 
which, respectively, couple to tp 1 -^ 2,3 for the electron mass and masses of atmospheric neutrinos 
and to ip 2 ' 3 il> 2 ' 3 for the fi- and r-masses and one-loop radiative neutrino masses relevant to solar 
neutrinos. This mechanism is realized by utilizing an additional residual discrete symmetry supplied 
by explicitly broken L', which guarantees the absence of tree-level neutrino mass terms of the 
ip 2 ' 3 ip 2 ' 3 -type. Pure rotation effects due to the diagonalization of neutrino and charged-lepton mass 
matrices are estimated to yield Am,Q/Am 2 tm £ (m e /m M ) 3 ^ 2 = 0(1O -4 ) but the radiative effects 
supersede the rotation effects to yield Amg/Am^ ira = 0(1O -2 ) as the LMA solution. 

PACS: 12.60.-i, 13.15.+g, 14.60.Pq, 14.60.St 
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I. INTRODUCTION 

In the recent report on the detailed analysis of solar neutrino oscillations jl|-|3| done by the Super-Kamiokande 
collaboration || , it has been suggested that solutions with large mixing angles are favored while solutions with small 
mixing angles are disfavored at the 95% confidence level. Neutrino oscillations controlled by such large mixing angle 
have also been observed and confirmed for atmospheric neutrinos Both observed oscillations are characterized 
by the same property of neutrino mixings, namely, large neutrino mixings, which are theoretically consistent with 
bimaximal mixing scheme J^J?]]. The difference arises in their oscillation scales denoted by Arn^ tm for atmospheric 
neutrinos and Atoq for solar neutrinos, which are specified by Am 2 tm ~3x 10~ 3 eV 2 and At/Iq ;$ 10 -4 eV 2 , thereby, 
indicating a hierarchy of Am 2 tm » Am|,. 

Neutrino oscillations arise for massive neutrinos |@]. Since neutrino masses of ~ 5.5 x 10~ 2 eV are implied by 
j, models for neutrino oscillations must be equipped with mechanisms generating such tiny masses of neutrinos 
which include breaking of the lepton number conservation for Majorana neutrino masses. It has been long 
known |l2| that such lepton number breaking interactions are inevitably contained in the minimal SU(3)l x U(1)n 
model for electroweak interactions (called the 331 model) |l3| because a charged lepton {€} and its antiparticle (£ c ) 
together with a neutrino (V) are placed in a triplet of SU(3)l as ij)={vi„ £l, ^1) T ■ ^ s a result, the model provides a 
charged-lepton mass term included in tpij), which simultaneously contains vl^l as a neutrino mass term. Namely, the 
need for Majorana mass term is linked to the existence of the massive charged leptons. Therefore, the 331 model is 
suitable for incorporating Majorana neutrino masses |12,nj-n6l. 
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With this nice feature of the 331 model in mind, we introduce the bimaximal mixing scheme into the model to 
accommodate observed neutrino oscillations. One of the possible scenario for bimaximal mixing scheme is to employ 
a new U(1)l> symmetry based on L e — — L T (= L') |l7],pl|. Neutrino mixings are predicted to be maximal for solar 
neutrinos with Arrig=0 but arbitrary for atmospheric neutrinos with Am^ tm ^ 0. The hierarchy of Am^ tm 3> A?71q 
is trivially satisfied since Am|=0. One has to next explain how a nonvanishing Am %. is induced while preserving 
Am 2 . tm Attiq . An elaborated mechanism is to generate Anig as radiative effects []19|] , whose generic smallness 
explains the relative smallness of Atoq over Am^ im ■ 

In this article, we construct our interactions that induce phenomenologically consistent neutrino and charged-lepton 
mass matrix, M„ and Mi, with appropriate radiative effects for Atuq. Since the model describes both matrices in a 
unified form, masses for charged leptons and neutrinos are correlated to each other. The mass term of ipip transforms as 
3* and 6, which requires 3 and/or 6* as Higgs scalars. The coupling of the triplet Higgs scalar to iptp is antisymmetrized 
with respect to the 5 , C/(3)i-index. The spin and statistics, then, requires antisymmetrized flavor indices carried by 
■0, in turn, generating antisymmetrized Mi with eigenvalues corresponding to — m T , which contradicts with the 
observed hierarchy of m p <C m T . The triplet Higgs scalar alone is not consistent with the observed lepton mass 
spectrum. Therefore, the sextet scalar, which provides symmetrized mass terms, is the key ingredient of the 331 
model ||]. [] 

Because the sextet Higgs scalar contains both an 5C/(2)i-doublet scalar for charged leptons and an 5{/(2)i-triplet 
scalar for neutrinos, we naturally expect "degeneracy" in the structure of their mass matrices. The charged leptons 
should contain a diagonal mass term such as the r-mass, which forces M v to have a diagonal term. However, M v should 
not contain diagonal terms to be consistent with bimaximal mixing scheme based on L' . To avoid this "degeneracy", 
we use the conservation of the L' quantum number supplemented by an appropriate Z 4 parity as well as a specific 
breaking pattern of L' yielding a residual discrete symmetry. The collaboration of these symmetries allows neutrinos 
to develop no diagonal mass terms. As a result, atmospheric neutrino oscillations are controlled by those tree-level 
off-diagonal masses while solar neutrino oscillations are to be induced by one-loop effects. ^] 

In the following section, we first introduce particles and their interactions and discuss how the required vacuum 
alignment is realized. Section III deals with the diagonalization of M v , which reflects rotation effects due to the 
diagonalization of Mi. The bimaximal condition on the neutrino masses determines a pattern of the matrix elements 
for the electron mass. One-loop radiative mass of v e is also calculated. The discussions in Sec. IV include to examine 
A?tJq by utilizing the correlation between M v and Mg and to confirm that the partial "degeneracy" between M v 
and Mi suppresses the rotation effects to yield the LMA solution based on the radiative effects. The final section is 
devoted to summary and discussions. 

II. MODEL 

The model contains neutrinos and charged leptons in triplets of ip z a — (y l L , £ l L , £ c l) T (h a=l,2,3) ^ with L' = 1 for ip 1 
and L' = —1 for ip 2,3 , where the superscript c of £^ L stands for the charge conjugation, and two (anti-)sextet Higgs 
scalars S^°'+^ with L'=(0, 2) for (S^°\ S(+)). 

Since S 1 - ' 1 can couple to ip 1 ^ -1 (i=2,3), z^l — v\ mixing terms are 
generated by (0|SW n |0) ^ and are to be responsible for the bimaximal structure. On the other hand, m T (and also 
m^) are generated by (0|S'^ + ^ 23 |0) ^ via its coupling to tp l ipi 2,3). Since ip 1 ^ also supply couplings to v^l 

and v t l, if (0|S'^ + ' 11 |0) ^ is also induced, then it disturbs the nice feature of the bimaximal mixing scheme based 
on V . Therefore, in order for the present approach to be acceptable, the simultaneous alignment of (0|S'^ + ^ 11 |0) = 
and (0|S l( ' + - )23 |0) ^ should be dynamically allowed. It is demonstrated that this vacuum alignment of S^ 11 and 
g(+)23 j g mc ieed possible to occur as a result of the emergence of a residual discrete symmetry specific to U(1)l>- 

A. Particles 

The 331 model is specified by the Z7(l)Ar-charge. Let Af/2 be the U(1)n quantum number, then the hypercharge, 
Y, is given by Y = — V3A 8 + N and the electric charge Q em is given by Q em — (A 3 + Y)/2, where A a is the SU(3) 



1 The use of an 5*(7(3)L-singlet vector-like heavy "electron" without introducing a sextet scalar |2^] has been advocated to 
describe consistent mass spectra of neutrinos and charged leptons |l6| , 

2 Similar discussions based on L' and a discrete symmetry have been done in Ref. |24| . 

3 Throughout this paper, the roman letters of i, j, ■ ■ ■ and the Greek letters of a, /3, • • • , respectively, stand for the three families 
and three S'l7(3)L-indices. 
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generator with Tr(A Q A h ) = 25 ab {a,b= 1...8). Let us first summarize the whole particle content in our model, where 
the quantum numbers are specified by (SU(3)l, U(1)n\ U(1)l')' 



^ = K,4^l) T :(3,0;1), V 



i=2,3 



(3,0;-l), 



(1) 



for leptons, 



Ql ={di,-ui,Jl) :(3*,-l/3;0), Q| = K,4,J|) :(3,2/3;0), 



1.2,3 



(l,2/3;0), d^ 3 :(l,-l/3;0), j£ a : (1, -4/3; 0) , ■/£ : (1, 5/3; 0) 



1.2 



for quarks 



r,= (f]°,T]-,r]+) :(3,0;0), p 
for triplet Higgs scalars, and 



(p+,p°,p++) J :(3,1;0), X =(X-,X 
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for (anti-)sextet Higgs scalars. These Higgs scalars have the following vacuum expectation values (VEV's): 
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and quarks and leptons will acquire masses via these VEV's, where the orthogonal choice of the VEV's in Eq.(||) and 
the vanishing VEV of (0|5 ( ^ + - )11 |0) in Eq.(^) are to be ensured by appropriate Higgs interactions. 

The pure S , J7(3)i-anomaly is cancelled in a vectorial manner. The anomalies from triplets of the three families of 
leptons and of the three colors of the third family of quarks are cancelled by those from antitriplets of three colors of 
the first and second families of quarks jljj. Other anomalies including U(1)n are also cancelled. The scale of the 331 
model is set by v x , yielding SU(3)l x U(1)n —> SU(2)l x U(1)y and we expect the magnitude of v x is of order TeV. 



B. Interactions 



To realize phenomenologically viable interactions, we further impose a Z± parity on the model, where Z±=-\- for 
ijj 1 , r\ and £W;=— for -0 2 ' 3 and S^;=i for p and x and similarly for quarks. Hereafter, we restrict ourselves to the 
lcpton sector. The Yukawa interactions for leptons are, then, controlled by the following lagrangian: 

-C Y = Yl fuWf^S (0) + £ \hWf^S (+) + (h.c), (7) 

i=2,3 i,i=2,3 

where /'s denote the Yukawa couplings with = fj^. The possible term of the ipiprj-type is forbidden by Z4. The 

Higgs interactions are described by self-Hermitian terms composed of 4> a 4> c ^ (<^=*7) P, x) and S^°' + ^ a ^ S^/^, , which 
include the potential term of V rjpx : 

Vvpx = x vp\v x p\ 2 + \ PX \p x X \ 2 + *xv\x x f]\ 2 , (8) 
where A's represent coupling constants and (a x b) a = e Q/37 a/3& 7 and by the non-self-Hermitian terms in 

V = popS^x + \iVvS (0)c S^ c + \ 2 T 1 S^p c x c + (h.c), (9) 

where fi denotes a mass scale and Ai 2 stand for coupling constants. There is a leptonic number associated with tp, 
which can be taken to be —2 for S^ ' + ^ and for all others. Since the Lg conservation is broken both spontaneously 
by S^ 0,+ ' and explicitly by Vq, there is no harmful Nambu-Goldstone boson. 
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Although the lepton number cannot be globally assigned to S , [/(3)i J -multipulets, to argue Majorana masses for 
neutrinos, it is useful to define the lepton number L to classify interactions p5| in the similar way used for the 
electric charge defined in the standard SU(2)l x U(l)y model. Following the classification lately done in Ref. [ p6| , 
the L-number can be introduced in the form of 

L = 4A 8 /2V3 + £, (10) 

where C = 2 for J^' 2 ; = 4/3 for X ; = 2/3 for Q*' 2 ; =1/3 for V* 1 ' 2 ' 3 ; = for u^ 2 ' 3 and d 1 ^ 2 ' 3 ; = -2/3 for Q\ , r), p and 
g(o,+). _ _2 f or j3 ^ which yields L = (1, 1, —1) for -0 1 ' 2 * 3 as usual. The assignment for leptons and Higgs scalars is 
shown in TABLE |. Among our interactions constrained by L', only r]S^p c x c breaks the L conservation; therefore, 
it is this interaction that gives the nonvanishing A few comments are in order to support our choice of VEV's of 
v n>PtX and v^J: 

1. The nonvanishing VEV of (0|5(°) n |0) is triggered by the L-violating rjS^ p c x° for A 2 < since (0|r?i|0) + 0, 
(0|p 2 |0) + and (0|xa|0) + 0; 

2. The nonvanishing VEV of (0|5 (0)23 |0) is triggered by the ^-conserving 

pS^x for m < since (0|p2|0) 7^ and 
(0|xs|0) ^ and by the L-conserving r]r]S^ c S Wc for Ai < since (0|t?i|0) ^ 0; 

3. The orthogonal choice of VEV's of rj, p and X as m Eq.dq) is supported by V TIPX if all A's are negative. It is 
because V npx gets lowered if 77, p and % develop VEV's. So, one can choose VEV's such that (0|?7i|0) 7^ 0, 
(0|p 2 |0) ^ and (0|x 3 |0) ^ 0. 

Thus, these VEV's are dynamically aligned by our Higgs potential. On the other hand, we cannot argue whether 
(0|5^ + ^ 11,23 |0) acquire VEV's or not since there are no triggers such as an L-violating r]S^ + 'r] and an L-conserving 
pS^Xi which can, respectively, align (0|S I ^ + - )11 |0) and (0|S'^ + ^ 23 |0) to be nonvanishing. However, so far these in- 
teractions are forbidden by the L' conservation. To invoke such dynamical alignment, thus, requires L'-breaking 
interactions. 

To introduce L'-breaking effects calls for careful consideration because of the specific alignment of VEV's of S^ + ' , 
i.e. (0|i?W 1:1 |0) = and (0|S' + ) 23 |0) 7^ 0. Its breaking readily induces dangerous interactions that creates (0|S( + ) n |0) 
such as tadpole interactions for S'' + ^ 11 . Owing to the vacuum alignment of (0|S'' + ' )11 |0) = 0, L'-breaking interactions 
should conserve L. There is one candidate that can keep this specific alignment, which is an L-conserving T)r]S^ c S^ + ' c : 

V b = X^S^S^ + (h.c), (11) 

where Ab represents a breaking of the L' conservation. Since U(1)l' is broken by the unit of |Ql'|=4, there is still 
a residual conservation due to a Z4 symmetry of exp(i7rQ£</2) to be denoted by Z^i, which is shown in TABLE || 
together with L' and Z4. This Zl> symmetry can be used to constrain possible forms of L'-breaking interactions. 
The mass terms for p and r are created by (0 1 5't+) 23 1 0) 7^ 0, which is now supported by Vb for A;, < 0. Therefore, Vb 
serves as a trigger for (0|5^ + ^ 23 |0) 7^ 0. 

The VEV of (0 1 11 10) would be created by the following L-violating interactions of 

• Tr(S(°)tS(+)) since (0|S(°) n |0) 7^0, 

• 7]S^r] since (0|r?i|0) ^ 0, 

. V S^p c X c since (0| m |0) 7^ 0, (0|p 2 |0) ^ and (0| X3 |0> 7^ 0, 

• p X S^ c S^ c and p X S(+) c S(+) c since (0|#j|0) 7^ 0, (0|xs|0) 7^ and (0|S(°'+> 23 |0) ^ and 

• Tr(S(°)S(°)sM), Tr(S^S^SW) and Tr(S^S^S^) since (0|S( a +> 23 |0) ^ 0. 

Fortunately, all of these interactions are forbidden by the interplay between Z4 and Zl>. Tabulated in TABLE |l|are 
Z4 and Zl> for non-self-Hermitian Higgs interactions allowed by SU(3)l x U(1)n, where the reader can recognize how 
these dangerous interactions are excluded. At this moment, we reach the plausible situation, where the alignment of 

(0|S (+)11 |0) = 0, (0|S (+)23 |0) 7^ (12) 

is dynamically supported. 
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III. MASS MATRIX 



The resulting mass matrices for charged leptons and neutrinos take the form of 

8m l l2 ^ m i3 





Mi = ( 5m{ 2 m| a m| 3 | , M v = | <Jm? 2 5mg d «Jm£f | , (13) 

TO 23 TO 33 



where are expressed in terms of VEV's: 



Smi^fuvjp, SmZi^fuvP, mf j = f ij v< 1 , +) , (14) 



leading to the proportionality of 



Sm^ cx Sm e li: (15) 



and 5mlj d stand for one-loop radiative masses to be estimated in the next subsection. We simply assume bm l u -C mjj 
(i,j=2,3) so as to reproduce the hierarchy of m e -C m^ T by the mechanism discussed in Ref. [E7J. 

Since the magnitude of neutrino masses Sm^ is controlled by , let us explain how naturally small vi°^ is obtained. 
An estimate of 

V^> ~ -X 2 Vr,V p /v x (16) 

can be found in the analyses of Ref. Jl5| applied to our Higgs potential. It calls for A2 < as has been expected in 
the vacuum alignment. The smallness of is, then, ascribed to either | A2 1 -C 1 |2(| or v x ^> v, hp pS[ ], the latter 



of which is known as the type II seesaw mechanism J29(] . where we use the mechanism due to | A2 1 <C 1 by simply 
expecting that v x is of order TeV. One will find that (— A2) ~a fewxlO -7 provides neutrino masses about 5.5 x 10~ 2 
cV. This is in fact a fine-tuning of A2 to reproduce neutrino masses so as to meet the value of 5.5 x 1CP 2 . However, 
we have already known this tiny magnitude, which is give by m e /m t . The relative smallness of m e is linked to the 
underlying chiral symmetry recovered at the limit of m e — > 0. This feature dictates the so-called "naturalness" fl30| . 
Similarly, since the L conservation is respected in the limit of A2 — > 0, the process of taking | A2 1 *C 1 is also regarded 
to be natural although both cases require other new physics to explain why such tiny values are realized. 

Our program to discuss patterns of neutrino oscillations is based on the use of the proportionality of Sm^ to Srriu 
indicated by Eq. ( |l5] ) . Since to generate the electron mass and to reproduce Am 2 tm will, respectively, determine 5m\ i 

and Smli , which yield the estimation of /v^ . The magnitude of Vu derived in Eq. ( |l6| ) now determines that of 



A2 once appropriate values of VEV's are given and, then, £>m™^- (i, j=2,3) to be defined in Eqs.(23) and (p4| ) are 
calculable for this estimated A 2 . In the end, solar neutrino oscillations controlled by Sm™^ are found to be consistent 



with the LMA solution. 



A. Charged Lepton Masses 

We first discuss the diagonalization of Mp. The following unitary matrix Up is obtained under the approximation 
of 8m l u <C m\j (i, j=2,3) and transforms Mi into M dla9 =lljMi Z7^=diag.(— m ei to m , m T ): 

1 - (6U 2 + 5U$) /2, SU^ SU T 

U t =[ - (c a SU^ + s a SU T ) c a (1 - SU 2 /2) - sJ 2 V Ttl s a (l - 8U 2 /2) + c a S 2 V^ ) , (17) 
SaSUf, - c a 5U r -s a (1 - SU 2 /2) - c a S 2 V T ^ c a (l - SU 2 /2) - s a S 2 V^ T 



up to the second order of SU^. T with 

= ((c 2 - s 2 a ) Sm[ 2 - 2c Q s Q 5mi 3 ) /A_, SU T = (2c a s a Sm[ 2 + (c 2 - s 2 ) Sm 13 ) /A H 
5% = X l SU l SU J /(X+ - A_), A± = (m e 22 + m e 33 ± ^ (m l 22 - m| 3 ) 2 + 4m< 3 2 ) /2, 



c a = cos a 



y /(m t 33 -X-)/(X + -X-), s a = sma = y/ (X+ - m^) /(A+ - A_), (18) 
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where i,j=fi, t for <5 2 Vy, A M = A_ and A r = A + . The masses of e, [i and r are given by 

m e = ATO12 + Ami 3 , m Al = A_+Ami2, m T = A + + Ami 3 , (19) 

where 

Ami2 = ((4 - s a) Sm{ 2 - 2c a s a 8m l l3 ) 2 /A_, Ami 3 = (2c Q s Q W 12 + (c 2 - s%) 5m[ 3 ) 2 /A+. (20) 

To be consistent with the observed mass pattern of the charged lepton mass matrix, we should present a mechanism 
generating plausible mass textures such as the hierarchical pattern J27J and the democratic pattern However, 
in the present context of the model, we are not aimming at clarifying all the physical aspect of the lepton sector 
but, instead, at clarifying the correlation of neutrino mass matrix with charged lepton mass matrix. Here, we simply 
parameterize mfj by adopting two distinct ways to realize m^ <C m r : hierarchical and democratic mass textures and 
discuss their influences in neutrino mixings. The resulting parameterization is given by 

1. the hierarchical pattern of 777,22 23 ^ TO 33 an d 

777 e w Sm^/m^ + Sm{ 3 /m T , m (1 sa m 22 + (m 23 /m T ) , m T ^ m e 33 - (m 23 / m T ) , 
8U^~5m 12 /m^, 8U T w Sm^ 3 /m T , (21) 

leading to cos a ~ 1 with sin a ~ ™23/ TO 33 1) an d 

2. the democratic pattern of m 22 rs 77123 ~ m 33 an d 

777 e w 5m 13 I trip + ^m^ 2 2 /m T , w Am, m r w 2t77o — Am, 

(5[/ M w —Sm^/m^, 5U T w 5m 12 /m T , (22) 

where mo and Am specify m 22 = = m-o an( i TO 23 = rn o — Am (m 3> Am), leading to cos a — sin a — 1/ y2. 
Therefore, we have found the estimate of <5mf 2 13 that reproduces m e . 



B. Radiative Masses 



Neutrinos acquire radiative masses of 8m\1 d i j (i,j=2,3) generated by one-loop interactions corresponding to di- 
agrams depicted in FIG.|l| and FIG.||. It is convenient to shift the basis for charged-leptons from the original one 
for Me to the intermediate one, where n and r have diagonal masses of A_ and A+, respectively. The corrections 
due to 5m{ 2 13 are neglected because of Sm[ 2 13 <C m^ r . In this basis, the couplings of S^ 3 ' to the lepton triplets 
are described by Futp 1 ^ 1 (z=2,3) with ('J 2 , ^ 3 )=(c a ip 2 — s a ip 3 , s a ip 2 + c a ip 3 ) and Fu are the rotated /ii-couplings 
similarly defined as (-F12, -Fi3)=(c a /i2 — Sa/13, Sa/12 + c Q /i3), which appear in the figures. 

The result of the calculation of 5m r 1 1 d i , (i,j=2,3) is given by 



Sm r ^ d 



2A 2 /J Y] F \i m ^ 



i=2.3 



Til 2 2 2\2,7-/2 2 

Hjn t i , m s(0)+ , m (0) _ , m + )v + I{m t i , m s(0)+ , 



,m 2 )v 



1 OA A .(0) ttt2 TV 2 2 2 2 \ i r / 2 2 2 2 \ 

+2A 2 Aitj^ 'F u m t i y(m ti , m s(0)+ , m s(0) _ , m, )+ ) + I(m^ , m s(0) _ , m s(0)+ , m J( _ ) 
£m™ d = 2\ 2 \ b vf ) FuFijmp 7(m 2 ;,m 2 (+)+ ,m 2 ( o)-,m 2 + ) + /(m 2 ; , m 2 (+) _ , m 2 (0)+ , m 2 _ ) 

A 



v n v p v x , 



for i,j=2,3 and m e =0 in this basis, with mp = A_ ~ m M and m^3 



J(a, 6, c, d) 
J(a, 6, c) 



I J(a,b,d) — J{a,c,d) 

16-7T 2 6 — c 

a In a 



61n6 



m r as in Eq.(18), where 



ilnt 



(a — 6) (a — c) (6 — a) (b — c) (c — 6) (c — a) 
The Majorana mass (5m™ d for u e is induced by effective couplings corresponding to FIG.|l|: 



(23) 
(24) 

(25) 
(26) 



G 



(W^ifcM^Ty, (27) 

for (j)=p or \, and 

e^7 e a^'7' V ,i V) i /(( g(+) ) ^ /(5 (o )) c^ (ffhriappx-r) (28) 

with similar terms scrambled by possible permutations of the fields. Both couplings conserve L' . On the other hand, 
the same effective coupling as Eq.(^8|) but with the replacement of A2 — > A&, thus of the L'-breaking type, yields 
8m r i °- d (i, j=2,3) for v^. T — v^, T , which correspond to FIG.||. 

C. Neutrino Masses 

Now, let us turn to examining the neutrino sector. The neutrino mass matrix is also transformed into M™ eak by 
Ue of Eq.([l7j) to maintain diagonal weak currents and is calculated, up to the second order of 5U^^ T , to be: 

/ SM 1U M 12 M13 \ 
M weak = uj MvUe ~ Ml2 S M 22 SM 23 , (29) 

y M13 SM 23 5M 33 J 

where 

M 12 = N 12 (1 - 2SU 2 - 8U?/2) - N 13 (5 2 V T ^ + SU^U T ) , 

M13 = ^13 (1 - 25U 2 - 8U*/2) + N 12 (S 2 V^ T - 5U^U T ) , 
SMu = SmZf - 2 (N 12 SU^ + N l3 5U T ) , 
8M 22 = 2N 12 6U„ + c 2 Jm r 2 a 2 d - 2c a s a Sm r 2 a 3 d + s 2 Jm r 3 a 3 d , 
SM 23 = N 13 SU^ + N 12 SU T + c a s a (6m™ d - 6m r 3 f) + (c 2 a - s 2 a ) Sm r 2 f, 

SM 33 = 2N 13 6U T + s 2 Jm r 2 f + 2c a s a Sm r 2 f + c 2 Jm r 3 f (30) 

-^12 = c a 8rrii2 — s q (5to^ 3 , N± 3 = c a Sm^ 3 + SaSm^. (31) 

It is obvious to see that M™ eak satisfies Tr(M™ eafc )=Tr(Af l/ )=c5m™< i + 8m r 2 f + 5m r 3 f and Tr(M™ eak T M™ eak ) = 
Tr(Mj ' M v )—2{8mi 2 2 + 8m\ 3 ). Since |Mi 2 ,i3| S> |5My|, the diagonal neutrino masses given by (m„i, —m v2l m v3 ) 
are calculated to be: 

m v \ = m v + m$ + m v2 = m v - m£f + mffl, m„ 3 = m^, (32) 

up to the second order of SMij equivalent to the second order of SU^t, where 

m„ = ^M 2 2 + M 2 3 , m$ = (SMu + cos 2 d v 5M 22 + 2cosz9„ sin$ u 5M 23 + sin 2 ® V 8M 33 ) /2, 

M' 2 M 12 

m ( Ji = sin 2 $ V 8M 22 - 2 cos $ v sin tf v 5M 23 + cos 2 $JM 33 , m {2) = + (33) 

2m v m v 

with cos "&V = M\ 2 jm v and sin d v — M\ 3 /m v as the mixing angle for atmospheric neutrinos and 
M[ 2 = (—(5 A/11 + cos 2 d y 8M 22 + 2co$,-d v sin $ V 8M 23 + sin 2 $ V 8M 33 ) /2, 

M' 13 = ((cos 2 d v - sin 2 ■&„) 8M 23 + costf„ sin$ u (8M 33 - 8M 22 )) jyft. (34) 
The masses of m^f and M' 12 13 can be further converted into: 

(i) _ Sm\f + 8mlf (M 13 8U^ - M 12 SU T ) (M 12 N 13 - M 13 N 12 ) 



with 



(35) 



M' M Ml MM Ml I ( A/ 12^M + ^13^r) (M 12 N 12 + M 13 N 13 ) 

M 12 = NiidUn + N 13 6U T H , (36) 



M '^T2 



NisdUfj, - N 12 dU T - 2 = 



(37) 
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where other irrelevant terms involving (5m™ d are omitted. The mass of m^ 3 is obtained by noticing the relation of 
Tr(M™ eak )=5m{f + 5rn^ d + <5m™ rf that yields 

2m« + m$ = 5m\f + 5m™ d + Sm™ d , (38) 

from which we find that 

(1), rod o (Mi3^ M ~ (M 12 jV 13 - M 13 iVi 2 ) 

m„ 3 — OTO33 z — 2 . (oyj 

To enhance the bimaximal structure in M™ eak is achieved by setting I-M12I ~ |M 13 |, where M 12 13 can be approxi- 
mated to be: 

M12 1=1 N12 = c a 5mi2 — s a 5mi 3 , M13 ss N13 = c a Sm" 3 + s a Smi 2 , (40) 

because I^L^rl *C 1. The condition of |Mi 2 | ss IM13I thus, requires, for the hierarchical case with c a ss 1 and s Q ss 0, 

<5r< 2 « (41) 

and, for the democratic case with c a = s a = 1/ a/2, 

3> <5"T-i3 or Sm^ 2 <§; 6m^ 3 . (42) 

The essence of the present research lies in the fact that the bimaximal mixing scheme for neutrinos, in turn, determines 
how the electron mass is derived. The bimaximal condition of Eq.(|H]) combined with Eq.(pl|) yields 



5m[ 2 ss 5m{ 3 ss ^/m e m M , <5C/ p ss ^Jm e /m^, 6U T ss ^m e /m T , (43) 
while Eq.([42"|) combined with Eq.(p2|) yields 

5m e 12 ~ cos (j)^m^m^ , Sm\ 3 ~ sin 0^/m e m^, <5£/^ ~ — svacj)^ m e /m^, 5U T ~ cos 0\/ 'm e /m T , (44) 

and 

<W[ 2 ss eiy/rriem^, 5m{ 3 ss ^/m e m^ STJ^ ss -^Jm^rn^ 5U T ss £i^/m^m^/m T , (45) 
where tan 2 </> -C to t /to m . 

IV. NEUTRINO OSCILLATIONS 

The neutrino oscillation parameters of Am 2 tm and Atoq are expressed in terms of our neutrino masses as 

Am 2 tm = m 2 2 - m 2 3 ss m 2 ss M 2 2 + Af 2 3 ss <5r< 2 2 + <5m^ 3 2 , (46) 
Am| = m 2 ! - m 2 2 sa 4m„m< 1 1 ) = 2m„ (tfrn^f + Sm^ + <5m™ d - m$) , (47) 



leading to 



Am atm J rad \ Am atm j r , 



with 

Ami \ _ 2 (Smlf + 6m r 2 f) ( Am% \ _ (M 13 ^ - j^ag^Q (M 12 jV 13 - M 13 iVi 



(48) 



(49) 



The ratio (Amg/AmJ (m ) ra j arises from the radiative effects and (AitIq/ Am^ tm ) ro t receives the first order contribu- 
tions from the rotation due to the diagonalization of M v ^. 

We are now in a position to estimating Arnf^ / Am 2 atm to meet the LMA solution. Before going to its explicit 
evaluation, we list our assumptions on numerical values of the relevant mass parameters: 



8 



• v x » v weak (= (2V2G F )- 1 /2 =174 GeV ) for SU(3) L x U(1) N -» Stf(2) L x 

• v n = v p — vf* = Vg = v weak /2 as the simplest case to keep v 2 + u 2 + vf^ 2 + v^ 2 = v^ eak for weak boson 
masses, 

• m n ± = m p + = m s (o)± = m s(+) ± = tWafc and m x + = |/x | = u x , 

where w x , presumably of order TeV, is left unspecified and is to be measured in the unit of 10v weak . 

A. Atmospheric Neutrinos 

Let us first reproduce Am 2 tm ~3x 10~ 3 cV 2 , which requires that 8m\ 2 + Sm^ 2 ~3x 10~ 3 eV 2 from Eq.(fl6"|). 

for Sm^, then, leads to 

4 0) w (5.5, 1.8/ cos 0,7.5) x 1CT 9 ^ 0) , (50) 

where three terms to be placed in the parenthesis, respectively, denote the cases with (hierarchical masses with 
<5m^ 2 ~ fimi3, democratic masses with bm\ 2 3> ^ m i3: democratic masses with <5m^ 2 "C 8m\%). In the following, the 
same notation is used. Since is expressed in terms of A2 and VEV's as estimated in Eq.(|l6|), we find the constraint 
on A2 given by 

-A a « (1.1, 0.36/ cos0, 1.5) v x /10v weak x 10~ 7 . (51) 

Therefore, the magnitude of A2 is not arbitrary but fixed to be | A2 1 = 0(1O -7 ) by VEV's in order to obtain Am 2 tm ~ 
3xlO~ 3 eV 2 . 



The relation of Sr 



at' 



W _ 



/vf> (i=2,3) together with Eqs.(43- 



B. Solar Neutrinos 



Next, we calculate Awig, which depends on the sum oi5m\i d and 8rn^ d given by Eqs.(23) and ( p4| ) with i = j = 2. 
To find numerical values of (5m™^-, we use the following property of the integral I of Eq.(25): 

1(0, a, a, b) re -l/a&(= -7i), 1(0, a, a, a) ps -l/a 2 (= -J 2 ). (52) 



By using the assumptions on mass parameters stated above, we know that the second term in Eq.(|23|) is proportional 
to I\ and all other terms including 8m r i ° ,d of Eq.(|24|) are proportional to Ii. In Eq. (^3|) , the first term thus dominates 
the second term because of I2 S> I\ for ™ 2 _ S> m2 (a)+ and v x ^> v p . The other terms proportional to I2, the third and 
fourth terms together with Sm^ , are more suppressed by the factor of Ai.t,v 2 yeafc /w 2 compared with the first term 
since \Xi,b\ are at most of order unity. Therefore, the first term in Sm\1 d gives the most dominant contributions to 
Attt-q and those from all other terms can be safely neglected. The radiative mass of Sm\i d saturated by the mediating 
r lepton contribution turns out to be: 

5m\f « -73A 2 (l, cos 2 </.m T /2m M , 1/2) (v x /10v weak f eV, (53) 

where, in the democratic cases, the factor of m T cos 2 (f>/m p is due to 8m\ 2 2 (o< /12) — m e m T cos 2 4> in Eq.(Q) and the 
factor 1/2 arises from c 2 = s 2 = 1/2. The rotated couplings of Fn used in these calculations are approximated to be 
(hi, /12/V2, e/13/V^) (i = 2,3) for 9=-l for i=2 and 6=1 for i=3. 

By inserting the above result of Eq.(|53]) into (Am 2 ^/ Am 2 tm ) ra d shown in Eq.(j4|) with m v ~ 5.5 x 10 -2 eV, we 
reach the following estimate: 

(Aml/Am 2 atm ) rad » -2.7A 2 (1, 8.5 cos 2 0, 0.5) (v x /10v weak ) 3 x 10 3 , (54) 

where m T /m p ps 17 is used in Eq.([5^). This estimate together with Eq.([5l"|) for A2 yields 

(Am 2 Q /A m 2 atm ) rad « (2.9, 8.1 cos0, 1.6) K/10zw fc ) 4 x 10~ 4 , (55) 

where cose/) is bounded from below by tan 2 <f> <C m r /m M . By taking 
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v x /10v weak « (l.8 - 3.3, (1.4 - 2.5) 2.0 - 3.6) , (56) 

we obtain Atoq = (10 _5 — 10~ 4 ) eV 2 relevant to the LMA solution to the solar neutrino problem and 

(5m™ d = Am|/2m l/ «10- 4 -10- 3 eV. (57) 

Since v wea k=^7A GeV, v x roughly lies in the range of 3 — 6 TeV, which is the right order of v x (= 0(1 TeV)) 
as anticipated for the scale of SU(3)l x U(1)n — > SU(2)l x C/(l)y. For v x — 4 TeV as a typical value, namely, 
v x /10v wea k=2.3, we find that 

(Am%/Aml tm ) rad *i (0.8, 2.3 cos </), 0.6) x lO" 2 , (58) 
leading to Am| - (2.4, 6.9 cos <f>, 1.8) x 10~ 5 eV 2 for Am 2 tm ~3x 10~ 3 eV 2 . 



C. Diagonalization Effect 

Finally, we consider effects from (Amg/Am^ m ) ro( , which may disturb the behavior specific to the LMA solution 

indicated by (Am|/Am^ tm ) ra( i. Since (Arag/Am^ m ) ro i is controlled by m£,g as in Eg. j4S|) , it can be found that 

the leading terms of 0(SU^. T ) contained in turns out to vanish. The cancellation can be seen from the insertion 

of the leading terms of M12 13 ~ -^12,13 into the equation of M 12j V 13 — Af 13 A r 12 appearing in the numerator of in 

Eq.(|39"|) that turns out to vanish. The next nonvanishing contributions to m^j are given by the second order terms 

in M12-ZV13 — M13N12 that result in the third order contributions to m„ 3 , which is beyond our approximation because 
possible third order terms for (i, j=2,3), which have not been included, directly give the third order contributions 

as can be seen from the definition of mj^ in Eq. (|35|) . Therefore, we expect that (Amg/Am^ m ) roj becomes at most 
0(6U*), which is about (m e /m A1 ) 3 / 2 (= 3.4 x 10~ 4 ) as shown in Eqs.©-©, leading to (Am 2 Q / Am 2 atm ) rot < 10~ 4 . 
This estimate of (Am,Q/ Am 2 tm ) rot should be compared with Eq.(^8|) for the radiative case and in fact gives more 
suppressed contributions. It is, then, concluded that we can safely neglect (Aitiq/ Am 2 tm ) ro t and that the model can 
provide the LMA solution to the solar neutrino problem as indicated by (Am 2 ^/ Am 2 atm ) ra d- 

Similarly, the third neutrino v 3 will also acquire masses at most of order m„(5[/ 3 , which is about 10~ 5 eV. The 
radiative correction to m„3 arising from 8m 33 d is estimated to be about 

Sm r 3 f ~ -0.2A 2 A b (1, m T cos 2 0/2m M , 0.5) (v x /10v weak ) eV. (59) 

Since A 2 ~ 10~ 7 and Af, < 1, we find that 5m 33 < 0.5 x 10~ 7 eV for v x ~ 4 TeV. Thus, altogether we observe that the 
bound on m„3 is given by m v3 55 10 -5 eV. 

The cancellation of the diagonalization effects on the neutrino masses of rn^i „ 3 is a reasonable consequence because 
the (1,2) and (1,3) entries of M v are exactly proportional to those of (see Eqs.(|63])-(|65|)) so that the rotation effects 
due to SUfj,. r , which induce SM22 = 2m y cosd v 5U^, <5M 23 = m v {cos , & v 5U^ + shi'& u 6U T ), 8M 33 = 2m„ sini9„<5[/ r , are 
partially cancelled. However, the diagonalization effects manifest themselves in the matrix element U e3 that connects 
u e to v z . The CHOOZ and PALOVERDE experiments imply that \U e3 \ 2 < 0.015 - 0.05 ||. For the mass matrix 
M weak of Eq.(||), this element is calculated to be 

U e3 = — — [(cos 2 § v - sin 2 <& v ) SM 23 + costf„ sm<& v (SM 33 - SM 22 )] 

= {(c a 5m" 3 + s a 8m\ 2 ) SU^ — (c a 5m\ 2 — s a 5m" 3 ) 5U T ] , (60) 

TTly 

leading to SU^/y/2 w (1, — sin<£, — l) x /m e /2m ll with m^/rrir -C tan 2 cj> <C m T /m p for Eq.(E3) allowing \SUfi\ ^> \8U T \, 
which includes the typical value of tan0 ~ 1. Our model predicts that |i7 e 3| 2 ~ (l,sin cf>,l) x 0.0025, which are 
consistent with the current experimental bound. 
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V. SUMMARY AND DISCUSSIONS 



We have successfully demonstrated that the minimal 331 model with the global L e — L^ — L T and Z4 symmetry, which 
involves two (anti-)sextet Higgs scalars, S^°' + \ indeed explains solar neutrino oscillations based on the LMA solution. 
The LMA solution is explicitly obtained by the use of the plausible parameter values for masses and couplings. The 
model has the following properties, which do not depend on the parameterization, that 

1. the global L e — — L T symmetry is explicitly broken by rjr]S^ c S^ + ' c down to a discrete symmetry defined 
by exp(i7rQ^//2) that forbids the L-breaking interactions for S^ + \ the absence of which is the essence of our 
discussions to be consistent with the bimaximal neutrino mixing scheme based on L e — — L T : 

2. the L conservation is explicitly broken by \2r]S^ p c x°, which generates the VEV of Vv^ for neutrinos calculated 
to be ~ — \2V n v p /v x , where A2 < is the trigger for ^= 0; 

3. the mechanism generating the diagonal electron mass via Sm^ (i=2,3) is linked to the requirement realizing the 
bimaximal neutrino mixing via Sm^ owing to the proportionality of Sm^ oc 5m l u ; 

and yields numerical results, which of course depend on the parameterization, that 

1. the smallness of m v \ « m u2 ~ 5.5 x 10~ 2 eV for Am^ m is realized by |A 2 | estimated to be 0(1O~ 7 ); 

2. solar neutrino oscillations are induced by 

• the one-loop radiative effects, which yield 

Am 2 Q /Am 2 atm ~ 1(T 2 , (61) 

corresponding to the LMA solution provided SU(3)l x U(1)n — > SU(2)l x U(l)y occurs at 3 — 6 TeV, 

• the effects of the simultaneous diagonalization of the neutrino and charged-lepton mass matrices, which 
are expected to yield 

Am|/Am 2 tm < (m e /m M ) 3 / 2 ~ 10" 4 , (62) 

and which certainly yield the U e s matrix element estimated to be \U e s \ ~ \/m e /2m fJi — 0.05, 

where, in the present model, the dominance of one-loop radiative effects yields the LMA solution. It should be 
stressed that the choice of v x =0(l TeV) automatically leads to the LMA solution. 

The present radiative mechanism generates the Majorana mass of v e , Sm\i d , given by the effective coupling of 
(rfip 1 )ip 1 x{xS < ' + * > ) c - This coupling is induced by rjS^p c x c and pS^x as i n FIG.|] and is of the L'-conserving 
type. The less-dominant radiative masses are induced by ip l ^p J ' (S^) c (S^) c r]px for v e v e with i,j=l contributing to 
8m\1 d and for v^^v^^ with i,j=2,3 contributing to 8m r i °- d . The i/^^i/^^-terms involve the L'-breaking coupling of 

rjr]S^ c S^ c as in FIG]|. The splitting of v\ and V2 is provided by Sm\i d , which is about 10~ 4 — 10~ 3 eV, leading 
to the LMA solution. The mass of 2/3 is provided by 5m\% at most of order 10~ 7 eV but the diagonalization due to 
M„ t i may give m„3 10~ 5 eV. 

We have also presented plausible order-of-magnitude estimation of neutrino oscillations based on M v , which is 
correlated to Mf. 

1. for the hierarchical texture with 8m\ 2 ~ 8m\^ 

(* ^/m f ,m tl y/mjn^\ / ^m e m M ^m e m M \ 

y/m^ * * , Mi = ^/rrTem^ 5M AM , (63) 

y/mjn^ * * / \ V^^V" M J 

2. for the democratic texture with 6m" 2 3> Sfn^ 3 , 

* c^^rriemr s^^/fa^m^ \ / c^y/m e m T s^m e m^ 

M v = e„ I C0 % /m e m r * * ,M = c^^/m^m^ M M — AM I (64) 

s^y/merrifj, * * J \ s </V me "V M ~ ^ M M 

with C0 = cos</> and = sin</> constrained by tan 2 <ft <C m T /m^, 
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3. for the democratic texture with <5m^ 2 <C Sr 



'13! 



M v = e v \ Eiy/mjmj; * * , M £ = e^m e m M M M - AM ] (65) 



with £<<1, 

where e v =Vv /vf 1 ', M ^> AM, £M and the asterisks stand for the suppressed entries to generate Anig. The partial 
"degeneracy" between M v to Me suppresses the rotation effects on Am| due to the diagonalization of M„ and Me. 
As a result, the solar neutrino oscillations are, respectively, characterized by AuIq ~ (2.4, 6.9 cos</>, 1.8) x 10 -5 eV 2 for 
v x ~ 4 TeV. Of course, if our mechanism explaining neutrino oscillations works further, one should next consider some 
reasons for the emergence of the approximate equality of &m\ 2 ~ ^ m i3 i n the hierarchical case and of the hierarchy 
of either Sm^ ^> Sm^ or bm\ 2 ^ ^ TO i3 m the democratic case. 

The masses of Higgs scalars are simply assumed to be about v wea k for those related to the weak boson masses for 
the order-of-magnitude estimate for radiative neutrino masses, namely, m v ± ~ m p + ~ m s (o)± ~ m s (+)± ~ v W eak- On 
the other hand, the \ scalar that generates masses for the exotic quarks and gauge bosons is assumed to have m x + 
~ v x (^> v wea k) because (0|x|0) is related to the spontaneous breakdown of SU(3)l X t/(l)jv to SU(2)l X U(l)y We 
have estimated v x to lie in 3 — 6 TeV in order to reproduce the observed neutrino mass pattern. It is further known 
that the 331 model contains peculiar particles, i.e. doubly-charged bilepton gauge bosons |33| with masses of the 
order ev x /^/2, that couple to £ c l J[i£l. Since v x can be 3 — 6 TeV, the bilepton gauge bosons can be as light as and 
even lighter than 1 TeV. It is, therefore, quite conceivable that we will observe the bilepton gauge bosons with masses 
around 1 TeV, which provide striking phenomena with clean £~-£~ jets pi[ |. 

Now equipped with the plausible mechanism that generates phenomenologically consistent masses for neutrinos and 
charged leptons, the 331 model may be a physically interested model of leptons. It not only explains three families 
of quarks and leptons because of the specific anomaly-cancellation mechanism used together with the asymptotic 
free condition on SU(3) C but also involves Majorana mass terms demanded by charged-lepton mass terms. There 
remain dynamical questions about the smallness of the L-violating coupling of A2 and the fine tuning of the lepton 
mass parameters to reproduce m e -C <C m t , which require some new physics beyond SU(3)l x U(1)l- It is also 
required to have extensive analyses on physics of the quark and gauge boson sectors as well to fully clarify potential 
power of the 331 model. 
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Table Captions 

TABLE H V ', Z4, and Zv for leptons and Higgs scalars. 
TABLE 0. Zi and Z v for non-sclf-Hcrmitian Higgs interactions. 
Figure Captions 

FIG. 1. One-loop diagrams for v £ v e with t 1 = fx (i=2) and r (i=3), where = and is the rotated /^-coupling 

defined by F12 = cos a/12 — sin a/13 and -F13 = sin a/12 + cos a/13. 

FIG. 2. The same as FIG.[j] but for v^.rV^^ with = v^l (i,j=2) and v t l (i,j=3). 



TABLE I. L' , £, Zi and Z^/ for leptons and Higgs scalars, where the lepton number L is related to C as L = 4A 8 /2y / 3 + £. 
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and Z L i for non-self-Hermitian Hig 


gs interactions. 
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FIG. 1. One-loop diagrams for V cL V ei with t = fl (i=2) and t (z=3), where m ti =F ii v^ ) and F 1; is the 
rotated / 1; -coupling defined by F 12 = cos a / 12 - sin a f n and F 13 = sin a / 12 + cos af n . 
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FIG.2. The same as FIG. 1 but for v^ LtL with v' L J =V ML (i ,j=2) and V TL (zj=3) 
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